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Abstract: Because of its acute sensitivity to distance in the nanometer scale, Forster 
resonance energy transfer (FRET) has found a large variety of applications in many fields 
of chemistry, physics, and biology. One important issue regarding the correct usage of 
FRET is its dependence on the donor-acceptor relative orientation, expressed as the 
orientation factor k^. Different donor/acceptor conformations can lead to values in the 
0 < < 4 range. Because the characteristic distance for FRET, Rq, is proportional to (k^)^^^, 
uncertainties in the orientation factor are reflected in the quality of information that can be 
retrieved from a FRET experiment. In most cases, the average value of k corresponding to 
the dynamic isotropic limit (<k^> = 2/3) is used for computation of Rq and hence 
donor-acceptor distances and acceptor concentrations. However, this can lead to significant 
error in unfavorable cases. This issue is more critical in membrane systems, because of 
their intrinsically anisotropic nature and their reduced fiuidity in comparison to most 
common solvents. Here, a simple numerical simulation method for estimation of the 
probability density function of for membrane-embedded donor and acceptor 
fluorophores in the dynamic regime is presented. In the simplest form, the proposed 
procedure uses as input the most probable orientations of the donor and acceptor transition 
dipoles, obtained by experimental (including linear dichroism) or theoretical (such as 
molecular dynamics simulation) techniques. Optionally, information about the widths of 
the donor and/or acceptor angular distributions may be incorporated. The methodology is 
illustrated for special limiting cases and common membrane FRET pairs. 
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1. Introduction 

Fluorescence spectroscopy has been utilized as a major tool in biophysics in general, and membrane 
studies in particular, for decades. Several parameters can be measured in steady-state (quantum yield, 
anisotropy) or time-resolved (intensity and anisotropy decays) conditions, and it allows monitoring of 
photophysical processes that are sensitive to distance, concentration (Forster Resonance Energy 
Transfer, FRET), aggregation (FRET, static quenching) or diffusion (coUisional quenching). 
Fluorescence's versatility and sensitivity have earned it a place among the most useful biophysical 
techniques [1]. In particular, FRET is a most useful approach to the investigation of diverse problems 
in membrane biophysics, including membrane protein mapping, lateral heterogeneity (membrane 
domains), determination of the transverse location (depth) of fluorescent residues/labels inside the 
membrane, protein/lipid selectivity (preference of a specific lipid for the protein vicinity), and 
membrane protein oligomerization [2], both in spectroscopic studies and, more recently, under 
the microscope [3]. 

Additionally, quantitative applications of membrane fluorescence studies to the recovery of 
structural and partition information frequently require, as input, results obtained from other techniques. 
Some of these can be measured experimentally, such as area/lipid, lipid molar volume or bilayer 
thickness. However, one parameter for which there is no experimental technique suited to a definite 
measurement (though it was shown by Dale and coworkers that intervals containing it can be inferred 
from adequate fluorescence anisotropy measurents [4]) is the FRET orientation factor, k . Because the 

2 1/6 

characteristic distance for FRET, Rq, is proportional to (k ) , uncertainties in the orientation factor are 
reflected on the quality of information that can be retrieved from a FRET experiment. Most often, the 

2 2 

theoretical value for the so-called dynamic isofropic limit (<k > = 2/3) is used, but <k > uncertainty is 
still widely regarded as an inconvenience [5] that may be especially important in membranes, because 
of their infrinsic anisofropic nature and the restricted rotational mobility experienced by fluorophores 
incorporated inside the bilayer. The 2/3 value corresponds to the djoiamic limit, where rotation of both 
donor and acceptor is fast compared with the excited state lifetime. 

As described in detail below in Section 3.1., the instant value of k for a given donor-acceptor pair 

2 

(which falls in the 0 < k < 4 range) can be calculated from adequate molecular frame vectors. 
Although this definition of k is not suited to experimental measurement, it can be conveniently used 
in molecular simulation studies. In particular, from the instant position coordinates of a molecular 
dynamics (MD) simulation trajectory, calculation of the orientation factor for a given FRET 
donor-acceptor molecular pair is straightforward, and averaging both over pairs and over time is 
conveniently carried out. This has been used in a number of literature reports (e.g., [6-9]) including 
membrane-embedded fluorophores [10]. However, this approach is not a choice for most experimental 
researchers, as setup and validation of MD simulations for novel fluorophores (for which 
parameterization is seldom available) is not trivial. 
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Even though computational calculations of are scarce, determination of membrane-embedded 



fluorophore orientations using computational (for reviews of MD studies of fluorescent membrane 
probes, see [11,12]) or experimental techniques (traditionally linear dichroism or other spectroscopic 
methods [13,14], or, more recently, methods based on polarized total intemal reflection fluorescence 
microscopy or other microscopic techniques [15-17]) are more fi-equent. 

This manuscript describes a simple manner to use data on probe location and orientation or 
alignment in membranes, obtained by experimental or molecular simulation methods, to provide 
estimates for the k distribution in the djaiamic limit, from numerical simulation. For this purpose, a 
simple program was developed to generate random planar location distributions of donor and acceptor 
probes. Fluorophore orientation is then generated for each donor or acceptor taking into account 
information provided by the user. Subsequently, k is calculated for a large number of donor/acceptor 

2 2 

pairs, and hence the k distribution and average value are obtained. Average k values for common 
membrane FRET pairs were obtained and are shown. For other FRET applications, the reader may use 
the software, which is available as supporting information to this article, without difficulty. This 
simplified approach enables immediate estimation of distributions for quantitative FRET 
applications, without the need to resort to sophisticated theoretical or experimental techniques. 

2. Results and Discussion 

2.1. Special Theoretical Cases 

This manuscript describes a simple manner to use data on probe orientation or alignment in 
membranes to estimate the distribution of the k FRET factor, using a small program described in more 
detail in Section 3. In this subsection we use the program to verify some results known from the 
literature, as well as to investigate other hypothetical scenarios of interest. 

2.1.1. Isofropic Dipole Distribution 

As described in Section 3, the software assumes normal distributions of cosOd and cosOa (truncated 
to the range -1 < cos0 < 1), where 0d and 0a are the tilts of donor and acceptor (respectively) relative 
to the normal to the membrane plane. The angles corresponding to the disfribution maxima are read, as 
well as the standard deviations and ob. If very wide angular disfributions are used (a » 1) uniform 
azimuth angular distributions = sin0/2 are in practice obtained, as illustrated in Figure lA. 
Therefore, and because cp distribution is uniform in the [0, 2k\ interval, isofropic orientation of both 
donor and acceptor is easily simulated. 

Figure IB shows the resulting probability density of k compared with the analytical result for 
isofropic dipoles, which is (e.g., [5]): 
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Figure 1, (a) Comparison of the simulated donor (red) and acceptor (green) density 
probabilities of orientation for planar distribution with the uniform distribution, sinO/2 
(black); (b) Comparison of the simulated density function (blue) with the theoretical 
curve of Equation 1 (black). <k^> = 0.66. Inputs: coplanar distribution of donor and 
acceptor, Od,max = 6a,max = 90°, oa = cjb = 3; A^= 10^ donor-acceptor pairs. 
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From the density distribution, one can calculate the average value, 



k') = jKV(K')^^K' =2/3 



(2) 



As shown in Figure IB, there is excellent agreement between both simulated and analytical 

2 2 

distributions and average k value. The result <k > = 2/3, described in the literature [18] for planar 
probe distribution, is identical to that for isotropic orientation of dipoles distributed in a three 
dimensional system. The latter can also be obtained from the program, with a slight modification for 
the position of the acceptor in each pair (allowing a randomly allocated acceptor z value — varying e.g., 
between 0 and a maximal value specified by the user — in each donor-acceptor pair). Again, there is 
excellent agreement with the theoretical expectations (not shown). 



2.1.2. Identical or Identically-Distributed Fluorophores 

This situation is illustrated in the inset of Figure 2. In this case, the distribution and average value of 
K for planar systems depend solely on the angle Od between the donor transition dipole and the normal 
to the membrane plane. An analytical solution to this situation was described by Knoester and van 
Himbergen [19], who showed that 



where 



2 2 , 
S + S^ 

3 3 



3(cos'^, 



(3) 



(4) 



is an orientation order parameter, equal to the average value of the second-order Legendre polynomial 
PiiQ)- For this particular case, since no orientation heterogeneity is (as yet) assumed, <cosOd> = cosOd. 
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This result, which arguably is little known in the membrane fluorescence research community, is 
however very pertinent e.g., for FRET among identical fluorophores (energy migration, energy transfer 
or homo-FRET). 

Figure 2 shows that the variation of <k^> against Od, obtained using our numerical approach, agrees 
with the analytical result. Large relative variations are observed, and between local maxima for Od = 0 
(<K^> = (invariant) = 1), 90° (<k^> = 1.25) and 180° (<k^> = (invariant) = 1), there are minima 
(<K^> = 5/9) at Od = acos((5/9)"^) = 41.81° and 180° - 41.81° = 138.19°. There are four Od values for 
which an average of 2/3 is obtained. These values correspond to the angles in the [0, 180°] interval for 
which S = 0 ^ cosOd = ± (1/3)"^ (the so-called magic angle 54.74°, and 125.26°) and 
5* = 2/3 <^^> cosOd = ± (7/9)"^ (28.13°, 151.87°). However, despite <k^> being identical to the isotropic 
value, the distributions are quite different. As shown in Figure 3, the distribution for 0d = 54.74° 
shows peaks at ~ 0.12 and ~ 1 (absent in the isotropic distribution), while it vanishes abruptly for 
K > 25/9 = 2.778. This latter higher limiting value corresponds to the situation of coplanar donor 
dipole, acceptor dipole and separation vector, with the two dipoles facing each other. On the contrary, 
in the isotropic limit, values of k near 4 (corresponding to coUinear dipoles) are allowed and no cutoff 
is observed. 

Figure 2, Average value of the orientation factor, <k >, as a flinction of the angle Od, 
between the fluorophore transition and the normal dipole to the membrane plane, for 
two-dimensional FRET between donors and acceptors with identical orientation (see inset 
for illustration). Calculations with Od values spaced by 7.5° were carried out to obtain the 
points (the error bars length corresponds to two standard deviations often short simulations 
performed with = 20,000 donor/acceptor pairs). The solid line is the analytical solution 
obtained by Knoester and van Himbergen [19]. 
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Figure 3, Comparison of the simulated k density function for Od = Oa = 54.74° (blue) with 
the isotropic theoretical curve of Equation 1 (black). 
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To illustrate the utility of this curve for quick estimation of <k^>, we consider the recent MD 
study of l-palmitoyl,2-[12-amino]dodecanoyl-5«-glycero-3-phosphocholine (C12-NBD-PC) in a fluid 
bilayer of l,2-dipalmitoyl-OT-3-glycerophoshocholine (DPPC, T= 323 K). The angular distribution of 
the transition dipole has a maximum around -120-125° relative to the bilayer normal [20]. From 

2 2 

Figure 2, this would be consistent with <k > = 0.67-0.77. Direct calculation of k using the atomic 
coordinates of the MD trajectory leads accordingly to <k^> = 0.72 ± 0.04 for molecules in the 
same-leaflet [10]. 

2.2. Effect of Orientation Heterogeneity 

One important feature of our numerical procedure for estimation of k distributions of 
membrane-embedded fluorophores is that it allows the introduction of orientation heterogeneity, in the 
form of a normal distribution of cosOd and/or cosOa, as mentioned previously. This type of 
heterogeneity is expected in simple rotational models such as the "wo bb ling-in-cone" motions [21]. 
The standard deviations of the cosOd and cosOa distributions (od and oa, respectively) are a measure of 
the extent of freedom of dipole orientation. Therefore, they serve a similar purpose to that of Dale et al. 's 
depolarization factors <du^> and <dA> [4]: a and <(f> are different ways to express the orientation 
distributions sampled by donor and acceptor. Whereas Dale et al.'s form of taking this orientation 
redistribution into account using depolarization factors was justifled by their perceived possibility of 
estimation of <k > ranges from time-resolved fluorescence polarization measurements, a is a 
technique-independent, more direct metric of the dipole orientation heterogeneity. 

Intuitively, upon permitting sampling of different donor/acceptor relative orientations, 
K distributions and average values are necessarily affected. As discussed above, the limiting case of 
a » 1 leads to the expected isotropic case. Here we address the more physical possibility of moderate 
heterogeneity (a ~ 0.05-0.3) and its effect upon the orientation factor. Figures 4 and 5 illustrate the 
alterations produced on distributions and average values upon widening the angular distributions of 
donor and acceptor dipoles for two distinct situations: identical donor and acceptor distributions (e.g.. 
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homo-FRET) centered around 0 = 90° (Figure 4) and quite distinct donor and acceptor distributions 
centered around 60° and 120°, respectively (Figure 5). In these calculations, a larger number of pairs 
(at least 10^) should be considered, to allow a better sampling of the Od and Oa distributions. In both 
cases, it can be seen that whereas moderate orientation heterogeneity produces visible alterations in 
distribution, significant effects require a very large degree of heterogeneity. This is more apparent for 
the system where donor and acceptor have very distinct angular distributions (Figure 5), where, 
remarkably, very wide distributions of both probes affect k by a degree <10%. 

Figure 4. Effects of widening donor and acceptor dipole angular distributions (identical, 
centered around 90°, depicted in (a)) on the simulated density flinctions (b); Colors 
indicate standard deviations of cosO distributions as follows: black, a = 0; blue, a = 0.05; 
red, a = 0.10; green, a = 0.20; magenta, a = 0.30. The numbers given in the inset of (b) 

2 6 

correspond to the average k value. A^= 10 donor-acceptor pairs. 
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Figure 5. Effects of widening donor and acceptor dipole angular distributions (distinct, 
centered around 60° (donor, sohd curves) and 120° (acceptor, dotted curves), depicted in 
(a)) on the simulated density flinctions (b); See Figure 4 for color scheme. The numbers 

2 6 

given in the inset of (b) correspond to the average k value. = 1 0 donor-acceptor pairs. 
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2.3. Distance Effects 

In most applications of intramolecular FRET (involving a single acceptor fluorophore as quencher 
of each donor), distance measurements are carried out assuming a single fixed value of Rq, assuming 
a <K > value (usually 2/3) independent of donor-acceptor distance. However, this view has been 



distance in a labeled protein [6], possibly related to constraints imposed by the linker used to attach the 
donor probe. Naturally, in systems where each donor may transfer its excitation energy to multiple 
acceptors at varying distances, correlation between donor-acceptor distance and their relative 
orientation will necessarily occur if the dynamic isotropic limit is not met. Consider for example that 
donors and acceptors are located in separate parallel planes (different depths in the membrane), and 
that both types of fluorophores have narrow orientation distribution. The range and relative weight of 
possible values will certainly differ for a donor-acceptor pair in which the acceptor is located 
directly below or above the donor, relative to the case where their lateral separation is large. Therefore, 
the question is not whether there exists correlation, but to which extent may vary as a function of 
distance, and how to best take this into account. 

For the sake of illustration, we consider the same donor and acceptor orientations as in the previous 
subsection, but assuming fixed (no orientation heterogeneity) 6d = 6a = 90° and (9d, 6a) = (60°, 120°). 
We also consider two separate possibilities in each case: (i) coplanar donor and acceptor spatial 
distributions; (ii) donor and acceptor planes separated by an arbitrary distance ofh = 2.0 nm. In each of 
these scenarios, the size of the square systems / was varied up to 1000 nm. 10^ pairs were simulated for 
proper sampling of donor/acceptor distances. It is clear fi-om the figure that location of both donor and 
acceptor dipoles in the same plane always lead to identical distributions (Figure 6A,E) and average 
values (Figure 6B,F). At variance, for donors and acceptors placed in parallel planes, there are very 

2 2 

significant changes in both k distribution (Figure 6C,G) and <k > (Figure 6D,H), when the size of the 
simulated system (which constraints the maximal possible donor-acceptor distance) is «100 nm. 

In this case, definition of an average becomes a problem. The best way to deal with this situation 
is to take notice of how FRET occurs from a given donor (with fluorescence lifetime in absence of 
acceptor equal to xo) to a distribution of M acceptors (where the distance between the selected donor 
and the ith-acceptor is denoted by Ri). The decay rate k is given by 



challenged, as MD simulations revealed considerable correlation between k and donor-acceptor 



M 



k — T Q H~ — 



(3) 



where the different individual FRET rates are calculated according to 




(4) 




Finally, the Forster radius for each donor-acceptor pair is calculated using 
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Figure 6. k distributions (a,c,e,g) and average values (b,d,f,h) obtained for identical 
(centered around 90°; (a-d)) and distinct (centered around 60° (donor) and 120° 
(acceptor); (g-h)) fluorophore dipole distributions, both for coplanar donors and acceptors 
(a-b,e-f) and with parallel donor and acceptor planes (separated hy h = 2.0 nm, (c-d, 
g-h)), for varying simulated system size, /. Colors in the left panels indicate the value of / 
(black, 1 nm; blue, 2 nm; red, 4 nm; green, 8 nm; magenta, 100 nm). All lines in panels (a) 
and (e) are essentially identical. 
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where Oq is the quantum yield of the donor in the absence of acceptor molecules, n is the refractive 
index of the medium, and the integral term is the overlap between the normalized donor emission (/(/I)) 
and the acceptor absorption {e{X)) spectra. The constant value in the equation above assumes nm units 
for both A, and Rq, as well as M~'cm~' units for e(A). Combining these three equations, one concludes 
that the FRET rate is given by 

M ^2 

^r=CX-^ (6) 

where C is a constant term for all acceptors in the distribution. That is, the natural way to average k 
over different acceptors, with different distances to a given donor, is to use the inverse sixth power of 
the donor-acceptor distance as weight. In this way, acceptors located more closely contribute 
signiticantly more to the average (as they are also responsible for most of the quenching by FRET). 
For the purpose of this calculation, we adapted our program for averaging in this way for FRET 
between a plane of donors and another of acceptors, and FRET in a bilayer geometry (with possibility 
of FRET between a plane of donors and two planes of acceptors — one in each bilayer leaflet). We then 
proceeded to calculate <k^> for the two angular distributions presented above (Od = Oa = 90°; and 
(Od, Oa) = (60°, 120°)) assuming a bilayer geometry with h\ = 0.0 and /?2 = 2.0 nm as the 
donor-acceptor interplanar distances. For coplanar fluorophores, an exclusion distance of 7?e = 0.8 nm 
was also introduced, as donors and acceptors in the same plane cannot lie closer than a given minimum 
distance value. This would be similar to a typical experiment, with acceptor distribution in both bilayer 
leaflets. For a better sampling of all possible orientations for small Ri values, a value of / = 100 nm was 
used. Ten different simulations were carried out in each case, to assess whether proper convergence 
was obtained. For 0d = Oa = 90°, an average of 1.237 with standard deviation 0.010 was obtained, 
whereas for (6d, Oa) = (60°, 120°) the corresponding values were 0.771 and 0.009. The relatively small 
standard deviations indicate that sampling was adequate. The average values were similar (very 
slightly lower) to the <k > calculated for the planar systems (1.250 and 0.776, respectively), and the 
small differences are due to the contribution of acceptors lying in the opposite leaflet. Because of the 
transverse distance of 2 nm, neglecting this leaflet does not lead to considerable error. In any case, for 
the following section that applies this methodology to computation of <k^> for common experimental 
donor-acceptor pairs, whenever possible, information on donor-acceptor interplanar distances was 
taken into account. 

2.4. Application to Experimental FRET Pairs 

Table 1 summarizes relevant information regarding membrane probes that are commonly involved 
in FRET experiments (as donors, acceptors or both), Forster pair combinations formed by these 
probes, examples of quantitative experimental FRET studies using these pairs, and information 
obtained from literature simulation or experimental studies regarding fluorophore location and 
transition dipole orientation of the probes. The latter information was used to calculate an average 

value using the inverse sixth power of the donor-acceptor distance as weight, as described in the 
previous subsection. These averages are also shown in the Table. It can be seen that all these values 
but one (homo-FRET between Rhodamine B probes) lie within 15% or 0.1 of the isotropic limit of 2/3, 
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even though the probe orientation distributions are quite varied. The most probable explanation for this 
is that all orientation data retrieved from the literature refer to liquid disordered bilayers, in which 
probes have generally broad angular distributions. Applications to FRET in liquid ordered or gel 
phases, where conformational freedom is largely reduced, would potentially lead to <k^> flirther from 
the 2/3 value. This is illustrated in the ?-PnA/DPH pair. For t-FnA, contrary to other probes, information 
regarding orientation distribution is also available in the gel phase (5°, 0.03 (gel) [22]). Using these 
input values instead of the fluid phase parameters, one would obtain <k > = 0.53 ± 0.01 instead of 
0.58 ± 0.01 for both probes in the fluid phase. It should also be emphasized that the dipole orientation 
details shown in Table 1 refer to a specific probe compound labelled with the fluorophore under 
consideration, and other probes exist where the same fluorophore could be expected to show different 
transverse location and/or orientation distribution in the bilayer. For example, NBD data are based 
in a MD study of NBD-PC, whereas FRET studies also employ other NBD probes (NBD-PE, 
NBD-cholesterol). 

2.5. Effect of <K> on FRET Efficiency 



The most common experimental observables in FRET are the decay of donor in the presence of 
acceptor, iT>\{t), and the FRET efficency, E. The decay law /da(0 for bilayer geometries and uniform 
probe distribution is available in the literature [23] 



'z)^(0 = exp 



exp 



nR^cy 



V^o y 



• exp 



7tR:c 



1-exp 



y 



yj 



(7) 



In this equation, c is the number of acceptors per unit area, y is the incomplete gamma function, and 
the other symbols have the same meaning as in the previous equations. Although Equation 7 was 
originally derived for a plane of acceptors containing the donor {cis transfer), it is also valid if the 
donor molecule is separated from the acceptor plane by a distance Re. From the time-resolved donor 
emission, the FRET efficiency, which is defined (and often experimentally measured) by 

E = '^-IdaIId (8) 

(where /da and /d are the donor steady-state emission intensities in presence and absence of donor, 
respectively), can be calculated using 



= 1 - j z^^ {t)dt / j z^ {t)dt 



(9) 



In the latter equation, /d(0 ^ exp(-^/To) is the donor time -resolved fluorescence decay in the absence 
of acceptor. 

To appreciate the effect of <k^> on the FRET efficiency, we flrst deflne i?o as the value of the 



Forster radius calculated with k = 2/3, and then consider a form of Equation 7 in terms of 
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dimensionless parameters ^= t/io (reduced time), o = n Ro c (number of acceptors in a circle of radius 



equal to i?o ), a= Rq /i?e and p = (3/2)<k^>: 

ioA (^") = exp(- ^)expi - p'^^cr 



r 



.1/3 



l-exp(-yg«Vy 



/ol/3 2 

p a 



(10) 



This form of the decay law is suitable for calculation (by integration over ^, done using Equation 9 
after substitution of variables t and Q of universal curves E vs. a, for fixed values of a and p. 
Representative examples are shown in Figure 7. 

Figure 7. FRET efficiency as a function of dimensionless acceptor concentration 
a = 71 i?o ^c, for given values of a = i?o IRa (panels ((a-d) correspond to a = 0.5, 1.0, 2.0, 

and 5.0, respectively) and <k^> (in each panel, curves obtained with <k^> = 0.30, 0.50, 2/3, 
0.80 and 1.25 are show in blue, red, black, magenta and green, respectively). 





From Figure 7, it can be seen that variations in FRET efficiency upon changing the <k > value used 
to calculate i?o are most important when the latter is smaller than or comparable to the exclusion 
distance, such as in panels (a) and (b). For example, for a = Ro /Re = 0.50, the FRET efficiency 

obtained for <k^> = 1.25 is >3 times larger than that obtained for for <k^> = 0.30 over the whole a < 5 
range. However, it must be noted that larger values of a = Ro I Re are often met in membranes, as long 

as i?o has a reasonable (-2.5 nm or larger) value and there are acceptors in the same bilayer leaflet as 
the donors. In these cases, such as illustrated in panels (c) and (d), the effect of <k^> on E, though still 
clearly detectable, is smaller, especially in the 0.50 < <k^> < 0.80 range which is observed for most 
common FRET pairs (Table 1). 
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Table 1. Donor and acceptor location and orientation literature data for common FRET pairs in membranes, and <k^> calculated from 
this information. 



Donor 
fluorophore * 


Acceptor 
fluorophore " 


Quantitative 
experimental 
FRET studies 


Donor transverse 
location 


Donor orientation 


Acceptor transverse 
location 


Acceptor orientation 

(."max? 


Calculated 

^^^^ <• 


t-VnA 


DPH 


[24,25] 


0.92 nm [26] ' 


20°, 0.15*^ 


0.75 nm [27] *■ 


0, 0.35 [27] '' 


0.58 ±0.01 


DPH 


NBD 


[28-31] 


0.75 nm [27] ^ 


0, 0.35 [27] ^ 


1.35 nm [20,32,33] ''^ 


125°, 0.35 [20] ^ 


0.66 ±0.01 


NBD 


NBD 


[34,35] 


1.35 nm [20,32,33]'' 


125°, 0.35 [20] ^ 


1.35 nm [20,32,33] ''^ 




0.66 ±0.03 


NBD 


Carbocyanine 


[36] 






1.26 nm [37]'' 


77°, 0.25 [37] '' 


0.76 ±0.02 


NBD 


Rhodamine B 


[25,35,38-42] 






2.1 nm [43,44] 


90°, 0.275 [43,45] « 


0.56 ± 0.02 


BODIPY 


Rhodamine B 


[30,46] 


0.8 nm [47] ^ 


10°, 0.5 [47] f 






0.57 ±0.01 


BODIPY 


Carbocyanine 


[48] 






1.26nm[37]f 


77°, 0.25 [37] '' 


0.72 ± 0.02 


Rhodamine B 


Rhodamine B 


[25,35,41,49,50] 


2.1 nm [43,44] ^ 


90°, 0.275 [43,45] « 


2.1 nm [43,44]'' 


90°, 0.275 [43,45] « 


1.08 ±0.01 



^ ^PnA: frawi-parinaric acid; DPH: diphenylhexatriene; NBD: 7-nitrobenz-2-oxa-l,3-diazol-4-yl)amino; BODIPY: 4,4-difluoro-3a,4adiaza-s-indaccnc; Relative to the center of the bilayer; 
° The values given refer to the B,^.,^ and ct values (read by the simulation program) that produce angular distributions most consistent with the simulation or experimental data reported in the 
cited references; Average and standard deviations of average for ten independent simulations with two planes, R\ ^ weight, N = \(f,R^ = 0.8 nm and / = 100 nm; ° Inferred from 
experimental fluorescence quenching data; ' Inferred from simulation data; * Inferred from linear dichroism data. 
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3, Methods 



3.1. Definition of k for a Given Donor-Acceptor Pair 



K is given by [5] 



= {cos I// J- -^cosy/^ cosy/^Y 



(11) 



where *Ft is the angle between the transition moments of the donor and acceptor and ¥d and *Fa are 
the angles between the donor and acceptor transition moments and the vector uniting their centers, R 



(Figure 8). An equivalent expression, involving scalar products of r = Rl 
ua (in the direction of the donor and acceptor dipoles, respectively), is 



R 



K =\Ud-Ua -3{r -UoJir -Ua) 



and unit vectors u d and 



(12) 



Figure 8. Vectors and angles relevant to the orientation factor. Adapted from [10]. 





Acceptor 



3.2. Calculation of k Distributions and Averages from Numerical Simulation 

Two programs, written as macros in the Visual Basic Editor within a Microsoft Excel document, 
were developed for this article. The programs are available as supplemental files to this article, and 
may be edited by any user. Microsoft Excel Visual Basic was chosen as it may be easily used by most 
readers, and produces immediate numerical and graphic output. 

Program 1 calculates distributions and average values without weighting, for a single acceptor 
plane. The basic algorithm can be summarized as follows: 

1 . Read inputs: OD.max, 6A,max, od, oa, /, N, donor-acceptor interplanar distance h, and the number of 

2 2 

bins (equal length subintervals of 0 < k < 4) for k distribution. 

2. Main cycle - calculate for the i"^ donor-acceptor pair: 

2.1. Set donor and acceptor coordinates, (xoi, ym, -^Di) and (xai, Jai, -^aO, respectively. For the x and 
j-coordinates, uniformly distributed random numbers between 0 and L are allotted, zoi is taken 
as zero, whereas zai is taken as h. 
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2.2. Calculate donor-acceptor separating vector R , and its modulus. 

2.3. Set donor and acceptor dipole orientation coordinates, (Odi, cpoi) and (Oai, cpAi), respectively. For 
(pDi and (pAi, uniformly distributed random numbers between 0 and In are allotted. Boi and 0Ai 
(the angles relative to the bilayer normal, so that they fall between 0 and tt) are calculated by 
firstly computing their cosines, using Gaussian random numbers centered around cos6D,max and 
cos6A,max, with Standard deviation od and oa, respectively. Gaussian random numbers are 
obtained using the Box-MuUer transform method [51,52]. If either of the calculated cosines 
falls outside the [-1,1] interval, new cosBoi and cosBAi are calculated. Otherwise, 0Di and Bai are 
obtained in the [0, tt] interval by inversion. 

2.4. Calculate Cartesian coordinates of unit dipole vectors ud and ua , fi"om (Boi, (poi) and 

(Oai, CpAi). 

2.5. Calculate k using Equation 12. 

2.6. Update histogram of k distribution (in fact, the program also calculates Od and 0a distributions, 
for testing and/or comparison and matching with distributions from the literature; these are also 

2 2 2 

updated at this point), and partial sum of k values (for calculation of <k > = sum(K i)/A^ at the 
end of the cycle). 
2.7 Increment i and repeat until i > A^. 

3. Write outputs to the active worksheet: average value and distribution, 0d and 0a distributions. 

All calculations and distributions shown in Figures 1-6 were obtained with this program. 

Program 2 follows a similar algorithm for calculation of for the i"^ donor-acceptor pair. However, 
it computes an average k value using Ri as weight for each pair. It may be used for both one and two 
acceptor planes. The program reads 0D,max, 0A,max, ^d, cJa, U N, and, for both planes, the respective 
donor-acceptor interplanar (hi and A2) and exclusion (i?ei and Rei) distances. Two main cycles are 
carried out, one for each acceptor plane, with calculation of distribution in each case. For both 
cycles, calculation of for the i"^ donor-acceptor pair is identical to that in the tirst program, with an 
additional check: if Ri < i?ei (first cycle) or Ri < Rei (second cycle), then new donor and acceptor 
distance coordinates are allotted and a corresponding new Ri is calculated. Partial sums of K^/i?^ values 

2 2 6 6 

are now computed, to obtain <k > = sum(K i/Ri )/sum(l/i?i ). As described in Section 3, this weighting 
scheme reflects the dependence of the FRET interaction with distance, and was used to compute the 
values of the right column of Table 1 (see footnote for typical input parameters). 

4. Conclusions 

A simple method for calculation of the distribution and average value of the FRET orientation 
factor K^, for membrane-located donor and acceptor fluorophores, is presented. The method is 
implemented in two programs that are available with the present article. The programs use as inputs 
the donor and acceptor transverse locations, orientation and orientational heterogeneity. Analytical 
results previously obtained in limiting cases were verified, validating our approach. A number of 
donor/acceptor orientation distributions tested in this work produce k distributions and average values 
which are significantly different to that corresponding to the isotropic limit, characterized by 

2 2 

<K > = 2/3. Considerable dependence of <k > on system size / when the latter is «100 nm was 
verified, which presents a problem regarding the most appropriate <k > to consider for the purpose of 
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quantitative analysis of FRET data. It is proposed that calculating <k^> using the sixth inverse power 
of the donor-acceptor distance is the most natural way to average k in this situation, and estimates of 
<K > were obtained for common membrane FRET pairs using available fluorophore location and 
orientation data. It was shown that for most of these systems, the calculated <k > lies within ±15% of 
the 2/3 value, possibly because probes in fluid bilayers have very broad orientation distributions. In 
these cases, substantial effects on FRET efficiency are not expected, provided that the Forster radius is 
significantly larger than the donor-acceptor exclusion distance. For more ordered bilayers, narrow 
orientation distributions may produce <k > further form the 2/3 value. In these ordered systems, an 
additional problem is that fluorophore dipole orientation randomization may occur in a time scale of 
the order or slower than the donor excited state lifetime. In case reorientation proceeds much slower 
than transfer, the opposite regime (static limit) is operative. As discussed in detail elsewhere [4,5], in 
this limit, it is generally not possible to use k averages to calculate donor-acceptor distances in 
intramolecular FRET (the only exceptions occur for isotropic dipole distribution or when the FRET 
efficiency is vanishingly small). In the intermediate case of comparable rates of reorientation and 
transfer, FRET kinetics becomes very complex [19]. Because these situations of restricted mobility are 
not contemplated in our method, which calculates distributions and averages corresponding to the 
dynamic limit, care must be taken when addressing these systems. However, even in such cases, the 
ease of application of this method — as long as fluorophore location and orientation data (or, at worst, 
"educated guesses") are available (in the literature or upon measurement) for the pair at hand — allows 
calculation of <k > which will constitute an improvement over the widely used but very seldom 
justified <K > = 2/3 assumption. 
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